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ABSTRACT 

The global nonlinear time dependent evolution of accretion disk-corona systems in 
X-ray binary sources has been investigated to provide an understanding of the low 
frequency (~ 0.04 Hz) quasi-periodic oscillations (QPOs) observed recently in the 
Rapid Burster MXB 1730-335 and in some black hole candidates sources (Cyg X-1 
and GRO J0422-|-32). We consider a-viscosity models in which the viscous stress is 
proportional to the total pressure. In contrast to previous time dependent studies it 
is assumed that all mass accretion and angular momentum transport take place in an 
optically thick disk, but that a fraction of the gravitational energy that is released is 
dissipated in a corona. It is found that the coronal energy dissipation can effectively 
reduce the amplitudes of the mass flow variations generated from the thermal and 
viscous instabilities (in comparison to models without a corona). Provided that the 
disk is close to a marginally stable state mild oscillatory nonsteady behavior results. 
These oscillations are globally coherent in the unstable regions of the disk. 

A model for the high and low states of black hole candidate systems is also 
proposed. It is suggested that the low state, which is characterized by a hard 
X-ray spectrum, corresponds to a disk configuration in which the inner disk is in an 
advection dominated hot optically thin state whereas the high state corresponds to a 
configuration in which the inner disk is in an optically thick state surrounded by a 
corona. In this model, the mass accretion rate in the system is higher in the low state 
than in the high state. The hard X-ray spectrum of QPOs observed in the low state 
can be naturally explained by such a model. 

Subject headings: accretion, accretion disks — black hole physics — instabilities 
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1. INTRODUCTION 

The spectral behavior of galactic black hole candidates (BHCs) and active galactic nuclei 
(AGNs) contains valuable information about the underlying physics of the accretion process. It is 
now well accepted that the soft X-ray black body spectral component from BHCs and the UV 
black body radiation from AGNs is due to a cold optically thick accretion disk. On the other hand, 
the power-law spectra from these sources suggest that optically thin hot matter is also present. 
The quasi-periodic oscillation (QPO) phenomenon in low mass X-ray binaries (LMXBs) and BHCs 
also provides a probe of the accretion process in these systems. In particular, an understanding of 
the observed luminosity variations may place important constraints on the physics of accretion in 
these objects. The QPOs in LMXBs are classified as horizontal-branch and normal-branch QPOs 
with frequencies in the range of ~ 20 — 50 Hz and 5 — 10 Hz respectively. The models proposed 
for these QPOs involve neutron stars (for a review see van der Klis 1989). Recently, low frequency 
QPOs with characteristic frequencies of ~ 0.04 Hz have been observed from the Rapid Burster 
MXB 1730-335 (Lubin et al 1992), Cyg X-1 (Angelini & White 1992; Vikhlinin et al 1992a, 1994; 
Kouveliotou et al 1992a), and GRO 10422-^32 (Vikhlinin et al 1992b, Kouveliotou et al 1992b, 
Pietsch et al 1993). In the Rapid Burster, the QPOs are observed in the soft X-ray energy band, 
whereas in BHCs the QPOs are observed in the hard X-ray spectral state of these sources. In 
the hard state of BHCs the source spectrum is characterized by a single power law and strong 
fiickering activity (see Tanaka 1991). The mechanism responsible for the production of these 
low frequency QPOs may be different from those proposed for the other classes of QPOs since 
these oscillations have been observed from both neutron star and black hole candidate systems. 
Milsom, Chen, & Taam (1994) and Chen & Taam (1994) recently proposed a disk instability 
model in which the thermal and viscous modes of the accretion disk are mildly unstable based 
upon a modified viscosity prescription, r = —ao(H / R)^pg, where r is the viscous stress, pg is the 
gas pressure, H is the local scale height of the disk, R is the distance from the central object, 
and ciQ and n are constants. This model was motivated by the fact that accretion disks based 
on a viscous stress proportional to the total pressure are unstable (Lightman & Eardley 1974) 
producing large-amplitude, burst-like fiuctuations (e.g., Taam & Lin 1984; Lasota & Pelat 1991; 
Honma, Matsumoto, & Kato 1991), whereas accretion disks based on a viscous stress proportional 
to the gas pressure are always thermally and viscously stable. 

More recently, Chen (1995b) suggested that a disk-corona system is naturally weakly unstable 
due to the stabilization effect associated with coronal energy dissipation confirming earlier work 
of lonson & Kuperus (1984) and more recently by Svensson & Zdziarski (1994). The disk-corona 
configuration is one of the general scenarios proposed to account for the presence of optically thin 
hot matter in X-ray binary systems (see, e.g., Liang & Price 1977). The concept of a corona 
has also been invoked in the interpretation of the spectra from Seyfert galaxies (see, e.g., Haardt 
& Maraschi 1991). Consequently, theoretical studies of the structure of a disk-corona system 
have been undertaken by a number of groups including Nakamura & Osaki (1993), Kusunose & 
Mineshige (1994) and Svensson & Zdziarski (1994). The inclusion of a corona leads to a reduction 
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in the ratio of radiation pressure to gas pressure within the optically thick disk as the fraction of 
the energy dissipated in the corona is increased. Hence, the mass accretion rate threshold above 
which the disk is unstable to thermal and viscous perturbations is raised (see lonson & Kuperus 
1984). 

In this paper, we report on the results of the time-dependent evolution of accretion disk-corona 
systems and investigate their possible relevance to the low frequency QPOs observed in LMXBs. 
In the next section we outline the assumptions and approximations underlying the model 
calculations. The numerical results of the calculations are presented in §3. Finally, the hard 
spectral characteristics of the observed QPOs in BHCs are incorporated into the framework of our 
disk instability model to provide a model for the high and low states of black hole transient binary 
systems in the final section. 



2. FORMULATION 



We consider a disk-corona system which is axisymmetric and non self-gravitating. The disk 
is optically thick and geometrically thin so that it can be described by the vertically integrated 
equations. The overlying corona is hot and optically thin. Following Svensson & Zdziarski (1994) 
it is assumed that all the mass accretion and angular momentum transport take place in the 
cold optically thick disk. We neglect the radiation drag which may be important in transporting 
angular momentum in the corona as in models of normal branch oscillations (Fortner, Lamb, & 
Miller 1989) in the bright neutron star X-ray binaries. In the case of BHCs, the importance of 
such effects is less clear. Since we focus on the thermal- viscous instability and do not consider 
the inertial-acoustic instabilities which have a much shorter timescale and may exist possibly 
in the innermost regions of the disk (see, e.g., Matsumoto, Kato, & Honma 1989; Nowak & 
Wagoner 1992; Wallinder 1991; Chen & Taam 1993, 1995), a Keplerian disk is assumed. In these 
approximations, the surface density of the optically thick disk, S, at a given cylindrical radius, r, 
can be obtained by combining the conservation equations of mass and angular momentum. It is 
given by the standard time dependent mass diffusion equation (Lynden-Bell & Pringle 1974) as 



9S _ 3^ 
dt r dr 



(1) 

where t is the time, and z/ is the effective kinematic viscosity. In equation (1), the Newtonian 
Keplerian angular velocity, flk = \/GM/r^, is implicitly assumed, where M is the mass of the 
central compact object. Relativistic effects on the gravitational potential are important near the 
inner edge of the disk, but we anticipate that the long timescale (~ 25 s) associated with the low 
frequency QPOs dictates that the region of instability lies at sufficiently large radii compared to 
the innermost disk that the results we obtain will be insensitive to the neglect of these effects. In 
addition, the lack of a strong soft X-ray emission component during the phases when QPOs are 
observed in BHCs suggest that an optically thick disk is not present in the innermost regions. The 
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viscosity is described by the standard a-model prescription as (Shakura & Sunyaev 1973): 

2 

V = -ac,sH, 



where a is a constant, Cg = ^P/S is the local sound speed, and H = Cg/flk is the half-thickness 
of the disk. Here P is the vertically integrated total pressure of the disk. 

The energy conservation equation for the cold accretion disk is represented by the balance 
between the local viscous heating, the local radiative cooling, and the global heat transport (radial 
advection). It is expressed as 



p 4-3/3 

r. - 1 



dliiT dlnT\ , , /91nS 91nS 91ni? 

Vr—^ - - 1) + V. 



dt dr J V dt dr dt 



(3) 

r or 

where T is the mid-plane disk temperature, [i = Pg/P is the ratio of the gas pressure to the total 
pressure, Fs = 1 -|- (4 — 3/3)(7 — l)/[/3 -|- 12(7 ~ ^)iP ~ 1)]? 7 = 5/3 is the ratio of specific heats. 
Here, the radial velocity, Vj., is given as 



Sri/2 dr 
and 

M = 2TrrT.Vr. (5) 

Although eq(3) is cast in a slightly different form than in previous studies, it is equivalent to 
the energy equation used in Chen & Taam (1994). Here, heating from the corona has not been 
included (see Svensson & Zdziarski 1994). This is justified by the fact that the heating of the disk 
by the corona only takes place in the surface layers down to a few Thompson optical depths (see 
Ross & Fabian 1993). Furthermore, it is known that the effect of irradiation is unimportant in the 
deep interior of the disk unless the irradiated fiux is comparable to the product of the viscously 
generated fiux and the optical depth of the disk (Fyutyi & Sunyaev 1976; Taam & Meszaros 1987). 
We assume that a fraction r] of the total gravitational potential energy released is dissipated in the 
corona (Haardt & Maraschi 1991). Since the physics of the dissipation in the corona is not well 
understood, r] is left as a free parameter. Thus, within the disk, the viscous heating rate per unit 
area becomes 

F+ = (l-^).s(r^)'. (6) 

In the optically thick radiative diffusion approximation, the local radiative cooling rate in the 
vertical direction is 

F- = (7) 

3kS ' ^ ^ 

where k is the sum of electron scattering and free-free opacities. The radiative energy fiux in the 
radial direction is written as 

Fr = -2HF-'-^. (8) 
or 
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Finally, for optically thick accretion disks, the equation of state is given by 

P = + 2i/— , (8) 

where is the mean molecular weight assumed to be 0.617. 

The optically thick disk is assumed to be truncated at an inner radius to allow for the 
possible existence of an inner optically thin region. We adopt such a picture to be consistent with 
the spectral data during the QPO active phase (i.e., the presence of a hard X-ray spectrum and 
the lack of a significant strong soft X-ray component). Thus, we postulate that the oscillations 
which originate in an optically thick outer cold disk, propagate to an optically thin hot inner disk 
and produce the observed X-ray QPO variability. 

The time dependent equations (1) and (3) are solved to calculate S and T via an explicit 
method. The accretion disk is divided into 41 grid points distributed equally on a logarithmic 
scale ranging from to an outer boundary, rgut = 50rj„. The initial structure of the disk is given 
by the thermal equilibrium solutions (see Chen 1995b). The outer boundary of the disk is chosen 
to correspond to a fixed temperature and surface density to ensure that the fiow of matter into 
the computational domain occurs at a fixed rate. At the innermost radius of the disk we impose 
the condition that the gradients of the temperature and surface density are zero. 



3. NUMERICAL RESULTS 



The accretion disk model is characterized by the mass of the central object, M, the mass 
accretion rate, M, the viscosity parameter, a, the coronal dissipation fraction, r], and the inner 
radius of the cold disk. We adopt the mass as IAMq for neutron stars and 10 Mq for black 
holes. For convenience, the mass accretion rate is expressed in units of the critical limit defined as 
Mc = 47r CM/ (KesCe), where Kes = 0.34 is the electron scattering opacity and e is the efficiency for 
the conversion of rest mass energy into radiation taken to be 1/6 for a neutron star and 1/16 for a 
black hole. The radii are measured in units of the Schwarzschild radius defined as = 2GM/c'^. 

To delineate the phase space in M, r], and for which the accretion disk is marginally 
stable, we have carried out time dependent calculations for a range of these parameters. The 
results are illustrated in Figure 1 for a black hole of mass IOMq and a parameter equal to unity. 
For greater fractions of the gravitational potential energy dissipated in the corona, the disk is 
stabilized after a few oscillations (less then 4-5), but at lower values of r] than that indicated by 
linear stability analysis (see lonson & Kuperus 1984; Svensson & Zdziarski 1994; Chen 1995b), 
due to the affects of advective energy transport. We point out that oscillations can occur from an 
initial steady state model due to the fact that the initial disk structure is not an exact solution of 
the time dependent equations. In particular, to construct the initial disk the radial diffusion term 
in the energy equation (see eq. 3) is not included and the advection term is only approximately 
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considered (see Chen 1995b). These terms act as initial perturbations to the disk which lead to the 
transient temporal development and eventually to either stabilization or oscillation. From Fig. 1 it 
can also be seen that the critical value r]c, above which the disk is stable, is a function of the mass 
accretion rate. In particular, r]c increases with the mass accretion rate for a given inner radius of 
the disk since a greater energy dissipation in the corona is required to compensate for a larger 
unstable region in the disk at higher mass accretion rates (in the absence of the energy dissipation 
in the corona). Furthermore, it can be seen that r]c asymptotically approaches a limiting value for 
M ~ Mc with the limiting value approached at lower mass accretion rates for smaller inner disk 
radii. For example, r]c approaches a value of 0.59 to 0.88 at M ranging from 0.7 — 0.3Mc for inner 
disk radii ranging from TOvg to lOr^. It is also found that r]c is greater for disks of smaller radii 
at a given mass accretion rate. This follows from the fact that the unstable region is larger for 
smaller inner disk radii, and, hence, greater energy dissipation in the corona is required to place 
the disk in the stable gas pressure dominant region. 

The time dependent calculations reveal that the characteristics of the global oscillations 
resulting from the development of the thermal- viscous instability are sensitive to parameters which 
deviate from those delineated by each curve, especially for those at large radii. For example, 
the oscillations are of greater amplitude for lower r] since a larger fraction of the disk is involved 
in the outburst. In addition, the timescales of the global oscillations are generally longer for 
disks characterized by larger inner disk radii and smaller fractions of energy dissipated in the 
corona. For oscillation timescales in the observed range of the low frequency QPOs (~ 20 — 30 s), 
parameters must be chosen which lie extremely close to a given curve. This would not likely to be 
the case in general and parameters deviating from these curves are expected. We point out that it 
is possible that r] may be affected by the instability in the disk. For example, for greater instability 
the coronal dissipation may be increased, thereby increasing r] and reducing the tendency toward 
instability. This may force the disk to always lie near the marginally stable state. Although this 
is plausible, it is difficult to quantify and, instead, we consider the simplest case that r] is not a 
function of the strength of the instability. In this case, it is found that significant changes in the 
properties of the oscillations would occur at large radii. On the other hand, the sensitivity of the 
results to the parameters deviating from the curve is much weaker for small inner disk radii (10 
rg - 30 rg) where the dissipation of gravitational potential energy in the corona must be greater 
to produce characteristics similar to the observed QPOs from BHCs. In the following we describe 
the numerical results from studies at these smaller inner disk radii. 

In the first sequence (see Table 1), which we have denoted as the standard model, we have 
followed the evolution of the thermal instability in a disk surrounding a black hole of IOMq 
accreting at M = 0.70Mc. It is found that the disk undergoes globally coherent oscillations. This 
behavior is illustrated in Fig. 2 where the evolution of the local temperature and column density 
is shown for four different radii. For reference we also include the steady state S curves which 
are locally thermally and viscously stable for ^ > and unstable for ^ < 0. It can be seen 
that the evolutionary paths take the form of loops with larger deviations from the steady state 
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structure present at smaller disk radii. Because of the importance of advective energy transport 
during the time dependent evolution, the paths in these diagrams during the unstable stages do 
not follow trajectories at constant S (see also, Chen & Taam 1994), but form loops centered near 
the lower turning point of the S curve. Note also that the average temperature is slightly less 
than the initial steady state one, especially at small radii. This is due to the additional cooling 
terms included in the time-dependent energy equation, namely, the radial radiative diffusion and 
the global advection. For this sequence the disk remains steady for radii > llAr-g. The timescale 
of the resulting mass flow variations is calculated to be 28 s with an amplitude of ~ 30% (see 
Table f and Fig. 3a). 

To determine the sensitivity of the results to the mass accretion rate, M was changed to 0.65 
and 0.75 in sequences 2 and 3 respectively, with all other input parameters held fixed as in the 
standard sequence. That is, we choose parameters which place the model disk nearly parallel to 
the marginally stable curve characterized by an inner disk radius of 30rp. The disk is very weakly 
nonsteady as shown in Figure 3b for sequence 2, but the amplitude and period of the oscillation 
for sequence 3 has increased slightly (Fig. 3c). Hence, the frequency of the oscillation increases 
whereas the relative mass fiow rate fiuctuation amplitude decreases as the mass accretion rate is 
decreased. However, this relation holds only with other disk parameters held unchanged and only 
over a limited range in mass accretion rates for which the disk is nonsteady. The latter refiects 
the fact that for higher mass accretion rates (> for these cases) the disk becomes stable due 
to the effect of radial advection (for example, see the upper stable branch in Fig. 2) 

The nonstationary behavior of the accretion disk is most sensitive to the coronal dissipation 
parameter, r]. For example, if we choose parameters away from the vicinity of the marginally 
stable curve by decreasing r] in the standard sequence to 0.7, the region where the disk is locally 
unstable increases and, thus, the disk exhibits larger amplitudes and longer period oscillations 
(sequence 4, Fig. 4a). On the other hand, if r] is increased above ric, i.e. if r] lies to the right of the 
curve at a given radius as illustrated in Fig. 1, the disk is stable. 

The global instability of the disk also depends on the location of the inner edge of the cold 
disk. For example, by decreasing Tin from 30 Vg to 20 r^, with other disk parameters held fixed, 
the spatial extent of the unstable region of the disk increases and the mass fiow modulations tend 
to have a larger amplitude and lower frequency (sequence 5, Fig. 4b) with respect to the standard 
sequence. By considering the infiuences of r] and Tin together, similar oscillation amplitudes and 
timescales to the standard sequence can be obtained with, for example, r] = 0.72 and Tin = 35rp 
(sequence 6, Fig. 4c). 

Both these properties of the oscillations are determined by the parameters of the disk-corona 
system. To be consistent with the frequency and amplitude of QPOs from BHCs, our calculations 
suggest that, a ~ 1, 77 ~ 0.75, Tin ~ 30rp, and M ~ 0.7Mc for a lOM© black hole. H the black 
hole is less massive, then a shall decrease approximately proportionally. We should also point out 
that, the result depends on M less sensitively for large r] then for smaller 77. Thus a larger r] > 0.5 
is preferred to avoid excessive fine tuning. A small value of a leads to lower disk temperatures 
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and increases the thermal and viscous timescales. In this case a smaller value of r] is required to 
destablize the disk. Thus we estimate a cannot be less than ~ 0.5. 

To determine the dependence of the oscillations on the mass of the compact object we 
decreased M to 1.4 Mq to model the evolution of an accretion disk surrounding a neutron star. 
In order to produce oscillations at frequencies ~ 0.04 Hz as observed in the Rapid Burster, the 
viscosity parameter, a, must be chosen to be smaller than the black hole case since the absolute 
size of the inner disk is smaller. The thermal and viscous time scales (relative to the dynamical 
time scale) would be correspondingly increased. Since a smaller a results in a cooler disk making 
it more stable, smaller values of both r] and are required to promote global instabilities in the 
disk. Figure 5 shows three examples which are relevant to the case of Rapid Burster. In all three 
sequences (7, 8, and 9), the fraction of energy dissipated in the corona is assumed to be 0.5. It 
can be seen that variability with frequencies in the range of 39 mHz to 64 mHz can be produced 
provided that the inner edge of the optically thick disk lies close to the neutron star surface. 
Specifically, the mass fiow rate at = 8rg in sequence 7 varies on a timescale of about 16 s 
with an amplitude of 50% for a mass input rate of 0.3 Mc and a = 0.02. Variations of the input 
parameters lead to similar trends as described for the black hole case. 

4. DISCUSSION 

It has been shown that low amplitude mass fiow variations can be produced in an optically 
thick accretion disk in which the viscous stress is assumed to be proportional to the total pressure. 
An essential feature in the model is the recognition that if a coronal region exists above the cold 
disk the nonlinear development of the radiation pressure induced thermal- viscous instability can 
be moderated. We have found that the stabilization of the disk depends not only on the amount 
of energy dissipated in the corona, but also on the location of the inner radius of the optically 
thick structure. Stabilization associated with r] had been noted by Chen (1995b, see also lonson & 
Kuperus 1984 and Svensson & Zdziarski 1994) whereas the stabilization associated with has 
been shown in the present work. This latter tendency toward stabilization can be understood as 
resulting from the reduction in the width of the unstable region. For a sufficiently narrow unstable 
region the disk is stabilized by advective energy transport (e.g., Chen & Taam 1994). These 
independent trends can be combined to illustrate the tendency toward stabilization. For example, 
for an inner boundary chosen closer to the central object, a greater fraction of gravitational 
potential energy must be dissipated in the corona to stabilize the disk. The numerical calculations 
demonstrate that the mass fiow variations can be of small amplitude provided that the fraction of 
energy dissipated in the corona, r], is high. Specifically mass fiow modulations on a timescale of 
20 - 30 s can be produced for r] > 0.75 for an inner boundary of the optically thick disk of < 30rg. 
For T] < 0.75 the amplitude variations increase in size (> 30%). The disk exhibits nonsteady 
behavior for a range in the mass accretion rates for a given r] and rin- For example, for r] = 0.75 
and rin = ^Or^ the disk exhibits mass fiow modulations for ~ 0.65Mc — lAfc for a equal to unity. 
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QPO behavior is expected for a modulation of the mass input rate into the inner disk region. The 
variation of the QPO frequency with respect to mass accretion rate is nonmonotonic since the 
inner radius of the optically thick disk is likely to vary as well. For higher mass accretion rates it 
is likely that the inner disk radius increases and these variations are expected to lead to opposing 
tendencies. That is, for a fixed inner boundary and a given fraction of gravitational potential 
energy dissipated in the corona, the QPO frequency decreases with increasing mass accretion rate. 
On the other hand, the QPO frequency tends to increase for larger inner disk radii at fixed mass 
accretion rate (provided that the inner radius is not increased to the extent that the disk becomes 
steady). This suggests that the variations of the QPO frequency depend less sensitively on the 
mass accretion rate than in the case where the inner edge of the optically thick disk is fixed in 
radius. 

The quasi-periodic oscillations observed from black hole candidate transient systems are 
characterized by a hard spectrum (e.g. Vikhlinin et al. 1992; Kouvelioutou et al. 1992). This 
indicates that the optically thin regions in the disk participate in the phenomenon. Such a disk 
structure is expected for a ~ 1 since the inner regions of such disks would be become optically 
thin for M > O.lMc. This disk configuration is similar to the two temperature inner disk proposed 
by Shapiro, Lightman, & Eardley (1976). In the present context, however, we assume that the hot 
inner disk is stabilized by the advection of energy and entropy into the black hole as conceived 
by Abramowicz et al (1988, 1995), Kato, Honma, k Matsumoto (1988) Chen k Taam (1993), 
Narayan k Popham (1993), Narayan k Yi (1994, 1995a,b) and Chen (1995a, b). In contrast to 
the Shapiro et al (1976) picture, we also require a corona in the outer cold optically thick disk as 
well to moderate the thermal- viscous instability. Thus, in our model, the inner disk responds to 
the mass fiow modulations produced by the thermal- viscous instability in the outer region. As a 
result, the QPOs will be characterized by a hard spectrum. Our model is also consistent with the 
ROSAT observation of QPOs of J0422-F32 in soft X-rays by Pietsch et al. (1993) provided that 
f'in ^ 30rp where the soft X-ray component can be produced in the optically thick disk. 

The low frequency QPOs observed in the soft X-ray energy band from the Rapid Burster 
MXB 1730-335 likely result from the modulations of the disk accretion onto the neutron star 
surface. The lack of hard X-rays suggests that either the electron temperature characterizing the 
optically thin gas is lower than in BHCs and/or that the optically thick disk extends closer to the 
central object. The fact that the electron temperatures in the optically thin regions of the disk are 
lower for neutron star binaries is expected because of the additional cooling provided by the soft 
X-ray fiux emitted from the neutron star surface (Sunyaev et al. 1991). The interpretation with 
regard to the presence of an optically thick disk close to the neutron star is also consistent with 
the fact that lower viscosities, which enable the optically thick disk to lie closer to the neutron 
star surface, are required to explain the timescale of the QPOs for lower mass systems. For these 
viscosities, the presence of an optically thin region is attributed to the existence of a small neutron 
star magnetosphere which disrupts only the very innermost regions of an optically thick disk. 

The transition from an optically thick to an optically thin disk surrounding a black hole at 
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can, in principle, be determined from theoretical considerations. It is known that there exists 
a relation between the mass accretion rate and disk radius for which an advection dominated 
optically thin solution is the only possible solution (Chen et al 1995) for accretion rates above a 
certain limit. On the other hand, for mass accretion rates below the limit, three solutions are 
possible of which two are stable (optically thick and advection dominated optically thin) and 
one is thermally unstable (optically thin solution in which local cooling is dominant). At these 
lower mass accretion rates, we suggest that the optically thick region extends to the black hole. 
Furthermore, we propose that the innermost regions of the disk are characterized by an optically 
thick structure surrounded by an optically thin hot corona. In this case, the hard X-rays are 
produced by a combination of the synchrotron process involving the relativistic electrons in the 
corona and the magnetic field in the disk and by inverse Compton scattering of the soft photons 
emitted from the cool disk. In this state, the soft X-ray photons result from the reprocessing of the 
hard X-rays produced in the corona which impinge on the cold disk and from the energy dissipated 
in the cold disk itself. We identify this state as the high state of BHCs. This interpretation is 
consistent with the observations of BHCs in which the soft X-ray spectrum is well fit by multicolor 
black body temperatures in the disk with an inner radius of 3 (Tanaka 1991; Tanaka & Lewin 
1995) and with the fact that the X-ray spectrum exhibits short timescale variability primarily at 
higher energies (e.g., Ebisawa 1991; Ebisawa et al 1994). In this model, the low hard spectral 
state of BHCs, wherein a significant fraction (> 99%) of the energy generated in the optically thin 
disk is advected into the black hole (Chen 1995a; Narayan & Yi 1995b), occurs at higher mass 
accretion rates than the high soft spectral state. Hence, the bolometric luminosity characterizing 
the low state can be lower than the high state even though the mass accretion rate is higher. 

Future theoretical investigations will be directed toward a more detailed study of the 
relation between the mass accretion rate and the inner radius of the optically thick disk. Such 
investigations are necessary for determining the relation of the QPO frequency with respect to the 
mass accretion rate and for delineating the mass accretion rate regime separating the high and 
low states of the BHCs. 

This research was supported in part by NASA under grant NAGW-2526. RET thanks Dr. F. 
Haardt for discussions while he was a visitor in Goteborg. 
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Sequence M(Mq) M(Mc) a 



1 10.0 0.7 1.0 

2 10.0 0.65 1.0 

3 10.0 0.75 1.0 

4 10.0 0.7 1.0 

5 10.0 0.7 1.0 

6 10.0 0.7 1.0 

7 1.4 0.3 0.02 

8 1.4 0.18 0.01 

9 1.4 0.3 0.01 



f'in{f'g) 


7] 


/(Hz) 


AM/M( 


30 


0.75 


0.034 


0.3 


30 


0.75 


0.038 




30 


0.75 


0.032 


0.5 


30 


0.7 


0.024 


1.1 


20 


0.75 


0.029 


1.0 


35 


0.72 


0.032 


0.3 


8 


0.5 


0.064 


0.5 


5 


0.5 


0.054 


0.8 


7 


0.5 


0.039 


0.55 



Table 1: Model Sequences 
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Fig. 1. — The M(r]) relation where the cold disk is marginally stable for a specified inner radius 
with M = IOMq and a = 1. Here the solid circles, solid squares, solid asteriks, and solid triangles 
correspond to disks with inner radii of 70 r^, 50 r^, 30 r^, and lOr^ respectively. For parameters 
to the left of each curve the disk exhibits nonsteady mass fiow modulations for a given inner disk 
radius. For parameters to the right of each curve the disk is globally stable for a given inner disk 
radius. The locations of model sequences 1, 2, 3, and 4 are illustrated for convenience. 
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Fig. 2. — Evolutionary paths of the temperature and surface density of the standard sequence 
for four radii. The dashed curves denote the steady state S curve solutions and are shown for 
comparison. Note that the evolutionary paths are described by loops with larger excursions from 
the steady state curve for smaller radii. 
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(c) Mo=0.75M^, rj^ = 30rg, 7^ = 0.75, a = l 
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Fig. 3. — Time variations of the mass flow rate relative to the input rate, Mq, at the inner edge of 
the disk. In all three sequences, M = IOMq, a = 1, t] = 0.75, and = SOr^. The top (a), middle 
(b), and lower (c) panels correspond to sequences where M = 0.7, 0.65, and 0.75 Mc- 
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(a) Mo=0.7M^, r;„=30r t] = 0.7, a=l 
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Fig. 4. — Variation of the disk mass flow rate at the inner edge of the disk relative to the input 
rate as a function of time. In sequences 4, 5, and 6 the central object is chosen to be a black hole 
with mass M = IOMq, the mass input rate is M = 0.7Mc, and a = 1. The top (a), middle (b), 
and lower (c) panels correspond to sequences 4, 5, and 6 respectively. 
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(a) Mo=0.3M^, rj„ = 8r 77 = 0.5, a = 0.02 
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Fig. 5. — Temporal variability of the disk mass flow rate at the inner edge of the optically thick 
disk relative to the mass input rate. In all three sequences, the central object is chosen to be a 
neutron star with mass M = IAMq and the fraction of energy dissipated in the corona is taken 
to be ?7 = 0.5. The top (a), middle (b), and lower (c) panels correspond to sequences 7, 8, and 9 
respectively. 



